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ABSTRACT 

Estimates indicate there may be ^ 10"^ low luminosity X-ray pulsars 
(L ^ lO"^"* erg s~^) in the Galaxy undergoing "low-state" wind accretion in Be/X-ray 
binary systems, and ~ 10® — 10^ isolated neutron stars which may be accreting directly 
from the interstellar medium. Despite their low effective temperatures (kTe ^ 300 eV), 
low luminosity accreting neutron stars with magnetic fields B ~ (0.7 — 7) X 10^^ G 
could emit a substantial fraction (0.5 — 5%) of their total luminosity in a moderately 
broadened (E/AE ~ 2 — 4) cyclotron emission line which peaks in the energy range 
~ 5 — 20 keV. The bulk of the thermal emission from these stars will be in the extreme 
ultraviolet/soft X-ray regime, which is subject to strong interstellar absorption and 
would be difficult to distinguish from spectra of other types of objects. In sharp 
contrast, the nonthermal cyclotron component predicted here will not be strongly 
absorbed, and consequently it may be the only distinguishing signature for the bulk of 
these low luminosity sources. We propose a search for this cyclotron emission feature 
in long pointed observations of the newly discovered candidate isolated neutron star 
MS0317. 7-6477, and the Be/X-ray transient pulsar 4U0115-|-63 in its quiescent state. 
We note that an emission-like feature similar to the one we predict here has been 
reported in the energy spectrum of the unusual X-ray pulsar IE 2259-|-586. 

Subject headings: accretion, accretion disks — line:formation — radiative 
transfer — stars :neutron — Xrays:stars 

1. Introduction 

Most of the 30 known X-ray pulsars in the Galaxy are bright, with luminosities 
L ^ Le = 1.4 X 10"^^ Mi.4 i?g ^Acap ((Tx/cm) erg/sec, where Le is the effective 
Eddington limit for magnetic polar cap accretion onto a neutron star of mass 
M = I.4M1.4M0, radius R = 10^i?6 cm, polar cap area Acap km^, and (Tx/'^m is the 
ratio of the nonmagnetic to the effective magnetic photon-electron scattering cross 
sections. Nearly half of these pulsars undergo wind accretion from Be star companions, 
but are usually detected only during high luminosity transient outbursts. These bright 
sources may just be the "tip of the iceberg", however, of a much larger underlying 
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population of lower luminosity magnetic accretors. Estimates suggest there are ^ 10"^ 
similar Be/neutron star binaries in the Galaxy accreting at low luminosities L « Le 
(Rappaport & van den Heuvel 1982; Meurs & van den Heuvel 1989; King 1991). 
In addition, assuming a production rate from supernovae of (10 — 100 yrs)~^ there 
should be ~ 10^-10^ isolated neutron stars in the Galaxy, of which some fraction may 
be accreting directly from the interstellar medium (ISM) at detectable luminosities 
(Trevis k Colpi 1991; Nelson, Salpeter, k Wasserman 1991; Blaes k Madau 1993). 

Although easier to detect, the emergent spectra from the high luminosity sources 
is extremely difficult to model owing to the strong coupling between radiation and 
the accretion flow (e.g. Arons, Klein, and Lea 1987). As a result, the theoretical 
interpretation of the observed spectra is not straightforward. By contrast, because 
radiation pressure is unimportant for low luminosity sources they could offer a much 
cleaner laboratory for testing the basic physics of accretion in a superstrong magnetic 
field environment. Here we present calculations of the emergent spectra of such low 
luminosity X-ray pulsars. In sharp contrast to their brighter counterparts, we find that 
sources with low effective temperature and field strengths B ~ 10^^ G should exhibit a 
unique nonthermal signature: a prominent cyclotron emission line, containing 0.5 — 5% 
of the total accretion luminosity, superposed on the Wien tail of the underlying 
soft thermal emission. If this line can be resolved it would enable an unambiguous 
determination of the surface magnetic field strength, thereby providing important 
constraints on neutron star magnetic field evolution. 

One can imagine a number of possible reasons why low luminosity accretion onto 
strongly magnetized neutron stars may not be able to occur, the most important of 
which is the centrifugal barrier at the rapidly rotating magnetosphere - the so-called 
propeller effect (lUarionov and Sunyaev 1975). Using standard spherically symmetric 
accretion theory, mass fiow onto the neutron star may be inhibited when the star 
rotates faster than the local Kepler period at the magnetospheric boundary. 

Pp ^ Pcr^t = Q0L-^^'i4oM-'J'R-^^'s, (1) 

where = 10"^°//3o G cm~"^ is the neutron star magnetic moment, and 
L = 10"^'*i34 erg s~^ is the luminosity from the neutron star. Equation (1) 
is valid for a neutron star moving relative to the local medium at speeds 
f ^ 4 X 10"^ km M^^J Rq^"^ L^^J fi^Q^"^ ^ so that gravitational focusing of the accreting 
fiow is important. It is not clear whether isolated neutron stars could have been spun 
down to these periods (Harding k Leventhal 1992; Blaes k Madau 1993). Nonetheless, 
it could be that these barriers merely reduce the accretion efficiency rather than 
preventing it altogether. Here we focus on what might be seen from a low luminosity, 
high magnetic field accreting neutron star if such systems can exist and are detectable. 
Clearly, the observation of a rapidly rotating (Pp < Pcrit), low luminosity pulsar would 
necessitate some refinement of current ideas on entry of plasma into neutron star 
magnetospheres . 
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In section 2, we describe our model for magnetic accretion and radiative transport. 
In section 3, we present spectra based on Monte Carlo cyclotron transport simulations. 
In section 4, we discuss possible applications of our model to the Be/X-ray binary 
transients with known magnetic field strengths and to a candidate isolated magnetized 
neutron star. A preliminary version of these results was presented at the Maryland 
Conference on the Evolution of X-Ray Binaries (Nelson, Wang, Salpeter & Wasserman 
1994a). 



2. Magnetic Accretion And Radiative Transport 



The physics of accreting matter decelerating in magnetized neutron star atmospheres 
is described in detail in Nelson, Salpeter, & Wasserman (1993) (hereafter NSW93). In 
strong magnetic fields, B ~ lO^^G, atmospheric electrons move freely along the field, 
but are quantized into Landau levels of energy En = (ra + 1/2)Eb perpendicular to the 
field, where 

Eb = — = 11.6^12 keV. (2) 

THeC 

Since radiative deexcitation greatly dominates coUisional excitation for typical 
atmospheric conditions and since ksT << Eb, atmospheric electrons normally reside 
in the Landau ground state (n = 0), and therefore constitute a sort of one dimensional 
gas (see, e.g., Ventura 1973). Consequently, the manner in which infalling matter 
decelerates and converts the bulk of its kinetic energy into radiation is dramatically 
different from what occurs in non-magnetic, three-dimensional atmospheres. 

Accreting protons, channeled along magnetic field lines, free fall into the 
plasma atmosphere at the magnetic polar cap and decelerate by undergoing 
Coulomb collisions with atmospheric electrons. In nonmagnetic atmospheres, 
these protons would penetrate to a stopping optical depth (Thomson units) 
T, = mpvjf/{6melnAc) = S.lMl^RQ^ilO/lnAc), where vjj = OMcMI^^R~^^^ is the 
proton free fall velocity and In Ac ~ 10 is the Coulomb logarithm for a nonmagnetic 
atmospheric plasma. In this case, all but a fraction ~ 1/lnAc of the initial proton 
energy is lost to electron recoil motion transverse to the proton velocity. In magnetic 
atmospheres, where electron motion perpendicular to the magnetic field lines is 
constrained, collisions with center of mass energy nigv'^ /2 > Eb can result in ^ ra 
Landau excitations up to a maximum level 
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NSW93 have shown that in moderately strong fields, where n^ax is fairly large, the 
stopping power of the plasma is well approximated by the nonmagnetic result, but 
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with a reduced Coulomb logarithm, 

InAc ^ln2nmax- (4) 

In particular, accreting protons penetrate the atmosphere to a magnetic stopping 
depth tb ~ Ts(ln Ac/ ln2nfnax) = ^^M^ ^R^^ / ln2nfnax and lose all but a fraction 
~ l/ln(2nmax) of their initial energies in exciting atmospheric electrons to 
1 < ra < rimax Landau levels. The decelerating protons do not veer appreciably from 
their original direction of motion along the magnetic field until almost all of their infall 
energy has been deposited. This is completely different from the stopping mechanism 
in strong fields, where Umax ^ 1 (e.g. Kirk & Galloway 1982, Miller, Salpeter & 
Wasserman 1987). 

In moderate magnetic fields, B12 ~ 1, where rimax is large, most of the incident 
energy of a decelerating proton is first converted into Landau excitations of the electrons 
it encounters along its path. An excited electron in the rath Landau level then decays 
back to the ground state primarily by successive electric dipole radiative deexcitations, 
thus emitting ra fundamental cyclotron photons via the cascade ra ra— 1, ra— 1 ra — 2, 
... 1^0. Emission of higher harmonic photons is rare mainly for two reasons: (1) 
the initial Landau excitations are dominated by ^ 1 transitions, and (2) the excited 
electrons are at most only mildly relativistic — their classical gyration velocities 
Vc = {2nEB I frieYl-^ ^ '^ff (cL eqn. [3]) correspond to Lorentz factors 7 ^ 1.3. (The 
situation is quite different for transrelativistic to very relativistic electrons [7 ^ 2] in 
strong magnetic fields, where multipole transitions, and hence higher harmonics, are 
important; e.g., Harding and Preece 1987.) For the parameters used here, NSW93 
found that initially only about 5% of the accretion energy is converted to cyclotron 
harmonics higher than the fundamental. The few higher harmonic photons produced 
have a large resonant optical depth (~ 7000(re/100 eV)~^/^rx), and quickly spawn 
fundamental photons via (predominantly) dipole radiative decay of the ra > 1 Landau 
levels they excite: the branching ratio of electric quadrupole to electric dipole decay 
rates is L„^„_2 /L„^„_i = {12 1 b){EB I me){n — 1) = 0.05i?i2(ra — 1). We do not follow 
higher harmonic photons in our numerical work; we assume that in a more detailed 
treatment, we would find that all such photons convert into fundamental photons 
almost immediately after birth. 

We emphasize that the cyclotron photons resulting from the cascade described 
above are produced nonthermally: their initial energy, Eb, may be far larger than 
the atmospheric temperature. Ultimately, the bulk of the energy in these "primary" 
photons is deposited back in the atmosphere as a result of electron recoil in magnetic 
Compton scattering, and photon absorption; this energy emerges in the form of a low 
energy continuum. The cyclotron photons that manage to escape form a nonthermal 
high energy spectral feature. 

We have used a Monte Carlo code to compute the polarized radiative transfer 
of cyclotron photons through the magnetized plasma (Nelson, Wang, Salpeter, & 
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Wasserman 1994b; hereafter NWSW). Cyclotron photons are born in the atmosphere 
with the profile derived in Nelson, Salpeter & Wasserman (1993). The dominant 
transfer effects are (1) magnetic Compton scattering which results in both angle and 
frequency redistribution (Comptonization), (2) polarization mode switching, and (3) 
absorption via inverse magnetic Bremsstrahlung. Of fundamental importance in our 
simulations is the low temperature of the atmosphere compared with the cyclotron 



where Fq = Laccrl Le- This has several consequences. First, the atmosphere may 
be treated as a cold plasma and we use the cold plasma polarization modes in our 
treatment of polarized transfer (Gnedin & Pavlov 1974, Ventura 1979, Nagel & Ventura 
1983). Second, unlike the case in hot magnetic atmospheres, the Doppler core of 
the resonant scattering line is narrow, so that cyclotron photons can quickly escape 
the line core and avoid resonant absorption (Wasserman & Salpeter 1980; Wang, 
Wasserman & Salpeter 1988; NWSW). This greatly increases the probability that 
photons escape the atmosphere. Finally, despite significant energy loss to Compton 
recoil, the escaping cyclotron photons are much harder than the typical thermal 
photon, so that the emergent cyclotron emission line will be prominently displayed 
against the soft continuum. In the limit when mode switching is rare, the line profiles 
may be obtained semi-analytically (NWSW). 



There are four parameters in our model: M,R,B, and Fq. Without exception, we 
take Ml. 4 = 1 = Rq. We therefore construct a two-parameter family of isothermal 
models in B and Fq. 

The shape of the emergent photon number spectra (angle averaged and summed 
over polarization) is shown in Figure 1 for B = 10^^ G (panel a) and 5 X 10^^ G (panel 
b). For both panels, Fq = lO""*; the spectral shape does not depend strongly on Fq. 
The continuum is assumed to be a black body spectrum with an effective temperature 
determined from equation (4). To this spectrum we add the emergent cyclotron photon 
spectrum, with the proper weighting, computed from the Monte-Carlo transfer code. 
The vertical arrows in the figure show the location of the cyclotron energy, Eb, where 
line photons are born. The photon energy degradation due to electron recoil (i.e., 
Comptonization) is clearly evident; despite originating at energies ~ 8 — 80 keV (for 
0.7 ^ Bi2 ^ 7), the emergent line centers lie in the range ~ 5 — 20 keV. 

Because of the strong energy dependence of the magnetic scattering cross sections, 
cyclotron photons escape in Ngsc r^/^ scatters rather than the ~ scatters expected 
in nonmagnetic Compton scattering. Consequently, although broadened and degraded 
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3. Emergent Line Luminosity and Spectrum 
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Fig. 1. — The emergent photon number spectra angle averaged, summed over polarization, and 
normalized to unit area from a magnetized polar cap neutron star atmosphere with (a) B = 10^^ G, 
Fo = Laccr/LE = lO'^, and Lune = 0.02i„,,„ and (h) B = SxlQi^G, = W~\ Lu^e = OmLaccr- 
The arrows indicate the cyclotron energy, Eb, where line photons are born. 

from their original energy, the accretion-induced cyclotron photons emerge as a clear 
line feature with E / AE ~ 2-4 for 0.7 ^ B12 ^ 7 (NWSW). The spectrum thus consists 
of a soft thermal component plus a hard nonthermal cyclotron emission line superposed 
on the Wien tail of the underlying thermal continuum. This general form should 
persist even for harder (e.g., power law) continua as long as the bulk of atmospheric 
electrons have energies ^ Eb (NWSW). 

Figure 2 shows the fraction of accretion luminosity that escapes in the cyclotron 
line, Liine/ Laccr ^s a function of magnetic field strength. The dependence on Fq enters 
only weakly through the absorption cross section, Gabs ^ l/^e oc Fq At high values 
of B, few cyclotron photons are produced so the line fiux falls off rapidly. The rapid 
fall-off at low B arises from the transfer microphysics: The energy lost per scatter by 
photons to electron recoil along field lines is small compared to the thermal Doppler 
width, so that line photons become trapped inside the line core where the scattering 
cross section is large (c/. Wang, Wasserman & Salpeter 1988; Wasserman & Salpeter 
1980) and are destroyed by absorption before they can escape. However, a significant 
line luminosity (0.005 ^ Lune/ Laccr ^ 0.05) is found for a broad range of field strengths 
0.7 ^ B12 ^ 7. In this regime, the field is sufficiently strong to avoid line trapping, but 
not strong enough to quench the initial cyclotron photon production. 

The soft X-rays below ~ 1 keV will generally be strongly absorbed by the 
intervening interstellar gas but may still be visible by instruments onboard satellites 
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Fig. 2. — The fraction of the total accretion luminosity in the nonthermal cyclotron emission line 
as a function of magnetic field strength. The curves correspond to Fq = (solid)10~^, (dashed)10~'*, 
and (dotted) lO'^. 

such as ASCA (0.5-10 keV) and ROSAT (0.2-2.4 keV). For an assumed continuum 
spectral shape, the soft X-ray component can be used to obtain the line-of-sight 
column density and the bolometric source luminosity. However, from the soft X-ray 
component alone, it will be difficult to distinguish the neutron star source spectra 
from spectra of other types of objects such as young stellar objects and active galactic 
nuclei. This is especially true of the ROSAT PSPC detector (0.2-2.4 keV) where the 
limited energy range above ~ 1 keV allows for a wide range of possible continuum fits. 
In sharp contrast, the nonthermal hard cyclotron emission line predicted here will not 
be strongly absorbed and thus serves as a unique identifying signature for strongly 
magnetized low luminosity accreting neutron stars. Note that while it is possible for 
the ASCA detectors to observe the nonthermal component (if B12 ~ 1), according to 
our calculations the ROSAT detectors cannot see this feature for any field strength. 



4. 



Discussion 



It may be possible to observe cyclotron emission line features in the spectra of 
highly magnetized X-ray transients during their quiescent phases, provided that some 
mass accretes onto these objects in spite of the propeller effect or any other factors 
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that may inhibit low luminosity accretion. For example, consider the Be/X-ray 
transient pulsar 4U0115+63 which has a rotation period Pp = 3.6 s, is relatively 
nearby at ~ 3.5 kpc (White, Swank and Holt 1983), and has a magnetic field strength 
of Bi2 ~ 1 as determined from cyclotron absorption features seen in its spectrum 
during transient outbursts (Makishima et al. 1990; Nagase et al. 1991) when 
Lfransient ~ 10"^'' erg s~^. If the propeller effect were 100% efficient then accretion 
could only occur for L > Lp^op ~ 3 X 10"^^ erg s~^. Suppose that in quiescence the 
star would accrete from the Be star wind at a luminosity L = fLp^o-p (with / < 1) 
if accretion were not inhibited at all, but, because of the propeller effect, only a 
fraction e << 1 of this amount actually accretes onto the neutron star. In this case its 
bolometric luminosity is Laccr ~ 3e/ X 10"^^ ergs~^, and, according to our calculations, 
the cyclotron line luminosity will be of order 50 times smaller (c/. Figure 2). A 40 
kilosecond integration with the Gas Imaging Spectrometer (GIS) onboard the X-ray 
satellite ASCA could obtain a Aa detection of the line feature in 4U0115-|-63 down 
to a fiux level of about 4 X 10~^'*erg cm~^ s~^ (between 5 and 10 keV). Thus, ASCA 
should be able to detect a cyclotron line feature if the quiescent accretion luminosity 
Laccr 3 X lO"^"^ erg s~^, or e/ ^ 10""^. Because of Compton recoil energy loss, this 
emission line will be centered at ~ 7 keV, well below the (fundamental) absorption-like 
feature seen at ~ 12 keV during transient outbursts. The Fe K emission line at 
about 6.4 keV, although nearly coincident with the cyclotron emission, should be far 
weaker and narrower. During outbursts, an unpulsed Fe K emission line has been seen 
centered at 6.4 keV with FWHM < 0.5 keV and i(Fe K)/ix(7 - 30 keV) ^ 9 X lO'^ 
(Nagase et al. 1991), and one might be concerned that this line could be confused 
with the nearby cyclotron emission feature during quiescence. However, we believe 
that the Fe K emission line should be far weaker than the cyclotron line. This Fe K 
line is ascribed to the reprocessing of ionizing photons ('^ 8 keV) by cool circumstellar 
wind material, leading to fiuorescent emission from Fe V - X. During quiescence, the 
thermal production of 8 keV X-rays is negligible, so the only potential source of 
ionizing photons is the cyclotron feature itself. Consequently, we expect the Fe K line 
to be ~ 100 times weaker than our predicted cyclotron emission line. 

The unidentified X-ray source MS0317. 7-6647 is a candidate for an isolated neutron 
star accreting directly from the ISM (Stocke et al. 1994). The position of this source 
coincides with that of a diffuse cloud seen by IRAS as an infrared cirrus (Wang and 
Yu 1994). If the source resides within this cloud, then it is at a distance ^ 100 pc. 
Its ROSAT PSPC spectrum (0.2 - 2.4 keV) is consistent with a 170 eV black body, a 
line of sight column density of Nh — 4 X 10^^ cm~^, and a total X-ray luminosity of 
2 X 10"^°((i/100 pc)^ergs~^ (Stocke et al. 1994). The inferred accretion patch area is 
only ~ lO""* of the total surface area of a neutron star, suggesting magnetic polar cap 
accretion. If this source is indeed an isolated neutron star and has a surface magnetic 
field of ~ 10^^, then a cyclotron emission line should be present between ~ 5-20 keV in 
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its hard X-ray spectrum. Detection of this line would strongly suggest that this object 
is an isolated, magnetized neutron star. 

Finally, we note that an emission-like feature at ~ 7 keV strikingly similar to the 
spectrum shown in Figure la has been seen in GINGA observations of the X-ray pulsar 
IE 2259-1-586 (Koyama et al. 1989). The statistical significance of this feature in the 
pulsed spectrum is marginal, so that further observations are needed to establish its 
reality. If it turns out to be real, then the feature could be due to Compton degraded 
cyclotron emission, as predicted here. If so, the surface magnetic field strength of this 
pulsar is ~ 1 X 10^^ G, and not the lower value ~ 5 X 10^^ G previously inferred. 

We thank 0. Blaes, Q. D. Wang, and A. Yoshida for useful discussions, and help in 
deriving the parameters for line detection by ASCA. This work was supported in part 
by NASA grants NAGW-666, NAGW-766, NSF grants AST 91-19475, AST 91-20599, 
AST 93-15375, and the NSERC of Canada. 
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